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Jamming transition in a highly dense granular system under vertical vibration 
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The dynamics of the jamming transition in a three-dimensional granular system under vertical 
vibration is studied using diffusing-wave spectroscopy. When the maximum acceleration of the ex- 
ternal vibration is large, the granular system behaves like a fluid, with the dynamic correlation 
function G(t) relaxing rapidly. As the acceleration of vibration approaches the gravitational accel- 
eration g, the relaxation of G(t) slows down dramatically, and eventually stops. Thus the system 
undergoes a phase transition and behaves like a solid. Near the transition point, we find that the 
structural relaxation shows a stretched exponential behavior. This behavior is analogous to the 
behavior of supercooled liquids close to the glass transition. 

PACS numbers: 45.70.Mg, 42.62.Fi, 64.70.Pf 



I. INTRODUCTION 

Granular materials are the particle systems in which 
the size of particle is large and the effect of thermal agita- 
tion is negligible 0. Recently there has been much inter- 
est in the physics of noncohesive granular materials lying 
on a vertically vibrating surface. When the vibration in- 
tensity is large, the granular systems show the properties 
of fluids, such as convection^ 0, heapingfl, trav- 
eling surface waves Jfj, pattern formation^], and size 
segregation^. When the vibration intensity is small, 
disordered granular materials become jammed, behav- 
ing like a system with infinite viscosity. After Liu and 
Nagel proposed an idea unifying the glass transition and 
the jamming behavior|lfjj|. the behavior of the jam ming 
transition has been studied extensivelv [llt Ha . Il3 . Hit Ha 
E3- IitL ni| . In this model, the inverse density of the 
system p -1 , temperature T, and stress a form the axes 
of a three-dimensional phase diagram, with the jammed 
state in the inner octant and the unjammed state outside. 
In athermal macroscopic systems like granular materials, 
thermal temperature does not play any important role. 
Instead, it is an effective temperature that relates the 
random motion of the particles. When jammed by low- 
ering the effective temperature, the system is caught in a 
small region of phase space with no possibility of escape. 
For thermal systems, if the molecules are bulky and of ir- 
regular shape, or if the liquid is cooled too rapidly for the 
crystalline structure to form, at low temperature it vitri- 
fies into a rigid phase that retains the disordered molec- 
ular arrangements of the liquid, creating a glass state. If 
the idea of the jamming phase diagram is correct, one 
can apply many ideas of the glass transition to explain 
the jamming behavior in athermal systems. 

In order to study the dynamics of granular systems, 
one might want to track the motion of the individual 
particles. However, due to the opacity of the granular 
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systems, most experimental work has been limited to 
the study of the external features of the granular flow 
or the motion of tracer particles. Recently noninvasive 
experimental techniques using magnetic resonance imag- 
ing (MRI), positron emission particle tracking (PEPT), 
and diffusing-wave spectroscopy (DWS) have overcome 
the problem of opacity in the three-dimensional granular 
system and have made observation of internal features of 
the flow possiblepa. Hfl 12(1 121| . Since the spatiotemporal 
scale of the particle fluctuations is much smaller than the 
resolution of MRI and PEPT, only DWS can provide a 
statistical description of a highly dense system of small 
particles with adequate resolution 22J. To briefly explain 
the DWS technique, photons are scattered consecutively 
by many particles in a highly dense medium. This dif- 
fusive nature of light in a strongly scattered medium re- 
sults in a scattered light intensity that fluctuates with 
time. Thus the temporal decay of the light intensity au- 
tocorrelation function is used to study very small relative 
motions of the scatterers in the medium. The DWS tech- 
nique has been successfully applied to study the motion 
of granular particles in a channel flow, a gas-fluidized 
bed, an avalanche flow, and a vibro-fluidized bed[2fl[2l|. 

In this paper, the fluidization and jamming process of a 
thick and highly dense three-dimensional vibro-fluidized 
granular bed is studied using DWS. Exploring the tem- 
perature axis of the jamming phase diagram [jjj, we com- 
pare our experimental result with theoretical concepts 
developed in the study of supercooled liquids close to the 
glass transition. Our result provides strong evidence of 
the analogy between the dynamics of granular materials 
and the behavior of supercooled liquids close to the glass 
transition. 



II. EXPERIMENTAL SETUP 

Figure 1 is a simple schematic diagram of the exper- 
imental setup. A 180 mm high, 90 mm wide (L x ), and 
10 mm thick (L y ) rectangular glass vessel is filled with 
glass beads of diameter 270 ± 20^im at the total depth 
of L z = 100mm. An electromagnetic shaker vibrates the 
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FIG. 1: Schematic diagram of experiment setup. 

vessel vertically with the form A sin (27r/t), where A is 
the amplitude and / is the frequency of the vibration. 
One can construct a dimensionless acceleration ampli- 
tude, whose maximum value is 

r = ^ 2 fA/g, (1) 

where g is the gravitational acceleration^. The vessel is 
much heavier than the total mass of the glass beads so 
that the collisions of the vessel with the granular medium 
do not disturb its vertical vibration. The air inside the 
vessel is evacuated below the pressure of O.lTorr, where 
the volumetric effect of the gas can be neglected [3 ■ The 
acceleration amplitude T and the vibration frequency / 
are monitored by the accelerometer which is positioned 
on the vessel. All the data are taken in a steady state 
condition where the system is well compacted [l5j. For 
the DWS measurements, an Ar laser beam with the wave- 
length of 488nm is directed on the wide side of the vessel, 
and the transmitted light intensity is detected by a photo- 
multiplier tube (PMT). The center of scattering volume 
is located at the position (x, z), where x is the horizontal 
distance from the central- vertical axis of the vessel and z 
is the vertical distance below the free surface. The inten- 
sity output / of the PMT is fed to a computer-controlled 
digital correlator (BI9000AT; Brookhaven Inst.). The in- 
tensity autocorrelation function, which is calculated by 
the digital correlator, is defined as 

g 2 {t) =< I(t)I* (0) > / < 1(0) > 2 = 1 + K G(t), (2) 

where t is a delay time, k is a factor determined by the 
optical geometry of the experiment |23| . and the bracket 
indicates a time average. In a highly dense granular 
medium, incident photons are multiply scattered by the 
glass beads. When the beads in the medium move, 
the intensity / measured at the detector fluctuates with 
time. Thus the intensity autocorrelation function has 
the information of the movement of the beads in the 
medium. When / is periodic, there are echoes in the 
correlation function at the integral multiples of the oscil- 
lation period|24|. In this experiment, the position of the 
incident laser beam is fixed and the scattering medium 
vibrates with the form of A sin (27T ft). Therefore when 
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FIG. 2: (a) The correlation functions measured using DWS at 
various dimensionless acceleration F in the case of the position 
(a; = Ornm, z = 50mm) and / = 50Hz. The unit of time in 
this figure is the vibration period T. (b) The height of the 
first echo in the correlation function at three different depths 
(z) in the case of x = 0mm and / = 50Hz. 

all the particles in the scattering volume exhibit perfectly 
reversible periodic motion, the correlation function will 
return to its initial value 1[=G(0)] at each multiple of the 
oscillation period. When only a fraction of the scatterers 
undergoes a reversible motion in the scattering volume, 
the height of the echoes will decay in time. Lastly, when 
all the scatterers exhibit an irreversible motion, there is 
no echo in the correlation function. Therefore the degree 
of fluidization can be characterized by the decay time of 
the height of the echoes in G(t). 

III. RESULTS AND DISCUSSION 

A. Fluidization point 

Figure 2(a) shows the correlation function at various 
dimensionless accelerations at the position (x = 0mm, 
z = 50mm) in the case of the vibration frequency of 
/ = 50Hz and the granular bed of L z — 100mm. When 
the maximum acceleration is smaller than the gravi- 
tational acceleration, T < 1, the correlation function 
oscillates with echoes at every multiple of the period 
T(= f )■ Here, the height of each echo does not change, 
that is G ec h (t) = 1. Since the photons meet the same 
scattering volume at every vibration period, G ec ho(t) = 1 
implies that there is no relative movement of the particles 
in the scattering volume. This suggests that the granular 
system is solidlike at T < 1, where the viscosity of the 
system appears to diverge. As the dimensionless accel- 
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FIG. 4: (a) The height of the echoes in the correlation func- 
tion up to the time of several hundred vibration periods. Here, 
the position is x — Omm, z = 50mm and / = 50Hz. (b) The 
same data in (a) with the scaled horizontal axis of (tjrf and 
the log scaled vertical axis. 



FIG. 3: The critical dimensionless acceleration F m at various 
position. Each data point is the averaged value from many 
independent measurements, (a) F m at various horizontal po- 
sition (x) in the case of z = 50mm and / = 50Hz. (b) Y m 
at various vertical positions (z) in the case of x = Omm and 
/ = 50Hz. (c) r m at various total depths L z in the case of 
the position (x = Omm, z = 50mm) and / = 50Hz. 



eration is increased above T = 1, however, the height of 
the echoes in the correlation function relax in time. At a 
certain critical dimensionless acceleration r,„, the height 
of the first echo disappears completely, [G(T) = 0]. For 
r > r m , G(t) decays very rapidly without any echoes. 
This implies that at T > T m , the particle positions be- 
come completely randomized by the external vibrations 
and the statistical description of the internal motion be- 
comes important. In this sense, T rn is the critical di- 
mensionless acceleration where the scattering volume be- 
comes fluidized. This fast decay in G(t) at T > T m con- 
tains information about the short-time dynamics of the 
granular particles^, [25|. Figure 2(b) shows the height 
of the first echo G(T) as a function of T at different verti- 
cal positions z. For a given position inside the vessel, the 
height of the first echo decreases monotonically with T. 
Thus the fraction of the particles which undergo an irre- 
versible motion increases with T. To determine T m , the 
data points are fitted with the function, tanh[a(r — r m )]. 

Figure 3(a) shows the critical dimensionless accelera- 
tion T m at various horizontal positions x with z=50mm 
and / = 50Hz. Note that T m is independent of the hori- 
zontal position. Figure 3(b) plots T m as a function of the 
vertical position z with x = Omm in the case of / = 50Hz 



and L z = 100mm. The critical value of the dimensionless 
acceleration increases linearly with the distance from the 
free surface, except near the free surface of the granular 
system. The solid line in this figure is a straight line with 
a slope of 0.0034. The result that T m is nearly 1 near the 
free surface of the granular system agrees with many pre- 
vious experimental studies 0, 0, 01 0- Figure 3(c) shows 
T m at various total depths L z in the case of the fixed po- 
sition (x = Omm, z = 50mm) . No dependence of T m on 
the total depth L z suggests that the relevant parameters 
in the fluidization process are just the distance from the 
free surface, z, and T m . According to this picture, when 
the maximum acceleration is larger than the gravitational 
acceleration (T > 1), there is a liquid-solid interface in 
the granular system. The upper part of the system is 
liquidlike and the lower part is solidlike, and the position 
of the interface depends on T. 

B. Relaxation near the transition point 

At 1 < T < r m , we observe a dramatic slowing down 
of the relaxation in the height of the echoes with decreas- 
ing r. Figure 4(a) shows the height of the echoes of the 
correlation function G ec ho (t) up to a delay time of several 
hundred vibration periods. The unit of time in this figure 
is the vibration period T. This shows that the dynam- 
ics of granular particles undergoes structural relaxation 
below r m . The data points are fitted with a stretched 
exponential decay function, 

G ech0 (t) = G exp{-(t/Tf}. (3) 
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FIG. 5: The characteristic relaxation time r of the height of 
the echoes from the same data in Fig. 4. 

Here, the characteristic relaxation time r depends on T. 
The fitting parameters of the stretched exponential decay 
function are G , r, and 13. During the fitting process, the 
value of G is fixed at G = 0.89, which is chosen from 
the data at T = 1.06 when the delay time approaches 
to 0. In the theory of the glass transition, the value 
of G depends on the scattering geometry and depends 
little on temperature and density [2(|. The fixed value 
of G is the major reason for the poor fitting at T = 
1.08. Figure 4(b) shows the same data with the scaled 
horizontal axis of (t/r) , where (3 = 0.50 ± 0.07. All the 
data collapse well on a straight line with a small value 
of the parameter (5 in Eq. ( 3) ( Here, (3=1 implies a 
simple exponential relaxation|23, Hg). This collapse is 
analogous to the behavior of supercooled liquids close to 
the glass transition. 

Figure 5 shows r at various T from the same data in 
Fig. 4. Here, the characteristic relaxation time r increases 
rapidly as T approaches a specific value. The solid line in 
this figure is a fit with the form of t ~ (T — r c )~ 7 , where 
T c = 1.048 ± 0.001 and 7 = 4.63 ± 0.06. The power law 
divergence of the characteristic relaxation time is similar 
with that of mode coupling theory (MCT). According to 
MCT, the characteristic relaxation time r at a temper- 



ature T in the glas sy relaxation process is of the form 
t ~ (T — X' c )~ 7 |26j. Here, T c is a dynamical crossover 
temperature from a liquidlike to a solidlike regime, and is 
located between the glass transition temperature T g and 
the melting transition temperature T m [Tj, |26( . From the 
analogy between the dynamics of the granular materi- 
als under vertical vibration and that of supercooled liq- 
uids close to the glass transition, one expects that the 
effective temperature of the granular system under verti- 
cal vibration is related to the dimensionless acceleration. 
Therefore we conjecture that r m , T c , and T = 1 in the 
dynamics of the granular materials under vertical vibra- 
tion correspond to T m , T c , and T g in the behavior of 
supercooled liquids close to the glass transition. 

IV. CONCLUSION 

The jamming process of a highly dense three- 
dimensional granular system under vertical vibration is 
experimentally studied using the DWS technique. At T > 
r,„, the granular system behaves like fluids. At T < 1, 
the system is arrested in the glassy state, where the vis- 
cosity of the system appears to diverge. At 1 < T < T m , 
the structural relaxation shows a stretched exponential 
behavior and the divergence of the characteristic relax- 
ation time r can be described by the power law function 
which is similar to MCT. This behavior provides the ev- 
idence of the analogy between the dynamics of granular 
materials and the behavior of supercooled liquids close 
to the glass transition. 
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